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Abstract 

Specialization to nectarivory is associated with radiations within different bird 
groups, including parrots. One of them, the Australasian lories, were shown to 
be unexpectedly species rich. Their shift to nectarivory may have created an 
ecological opportunity promoting species proliferation. Several morphological 
specializations of the feeding tract to nectarivory have been described for par- 
rots. However, they have never been assessed in a quantitative framework con- 
sidering phylogenetic nonindependence. Using a phylogenetic comparative 
approach with broad taxon sampling and 15 continuous characters of the diges- 
tive tract, we demonstrate that nectarivorous parrots differ in several traits from 
the remaining parrots. These trait-changes indicate phenotype-environment 
correlations and parallel evolution, and may reflect adaptations to feed effec- 
tively on nectar. Moreover, the diet shift was associated with significant trait 
shifts at the base of the radiation of the lories, as shown by an alternative statis- 
tical approach. Their diet shift might be considered as an evolutionary key 
innovation which promoted significant non-adaptive lineage diversification 
through aUopatric partitioning of the same new niche. The lack of increased 
rates of cladogenesis in other nectarivorous parrots indicates that evolutionary 
innovations need not be associated one-to-one with diversification events. 
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Introduction 

Although most flowering plants are pollinated by insects, 
a considerable number of tropical angiosperms are polli- 
nated by birds and bats specialized on nectarivorous diets 
(Bawa 1990; Sekercioglu 2006; Fleming and Muchhala 
2008). The associated ecological specializations resulted 
in several radiations of nectarivorous birds and bats in 
the tropics and subtropics (Fleming and Muchhala 
2008). Nectarivory has evolved convergently in several 
groups of birds, with the Neotropical hummingbirds 
(Trochilidae, 325-340 species), the Australasian honeyeat- 
ers (Meliphagidae, 182 species) and the sunbirds 
(Nectariniidae, 132 species) of Africa and Australasia, 
representing three major radiations of nectarivorous 
birds. Additionally, nectarivory can also be found in 



several groups of parrots (Psittaciformes). Parrots represent 
one of the most species-rich clades of birds (Jetz et al. 
2012). WhUe they feed mainly on seeds and fruits, the 
chiefly Australasian Loriinae (lories) are specialized on a 
nectarivorous diet (Collar 1997; Rowley 1998). The lories 
consist of 53 species (Collar 1997) and are considered 
generalized flower visitors with eucalypts being a particu- 
larly important nectar source (Fleming and Muchhala 
2008). Besides the lories, the swift parrot Lathamus discolor 
of Australia, the genus Loriculus of Australasia and Indo- 
Malaysia as well as the genus Brotogeris of the Neotropics 
are all supposed to depend on nectar as food (Homberger 
1980; Giintert 1981; Forshaw 1989; Collar 1997). Their 
specialization to nectarivory has evolved in convergence 
to that of the lories (cf Wright et al. 2008; Schweizer et al. 
2010, 2011). 
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Nectar is a liquid food source rich in sugars, which 
account for almost 100% of its dry weight (e.g., Liittge 
1976; Gartrell 2000). However, it contains only a small 
amount of amino acids, too low to satisfy the nitrogen 
requirements of a bird (Martinez del Rio 1994). Therefore, 
nectarivorous birds have to rely on other nitrogen sources 
like insects or poUen (Richardson and WooUer 1990; Brice 
1992; van Tets and Nicolson 2000; Nicolson and Fleming 
2003). Several morphological and physiological specializa- 
tions to nectarivory have been described in the various nec- 
tarivorous bird groups (Schlamowitz et al. 1976; Brown 
et al. 1978; Richardson and WooUer 1986; Casotti and 
Richardson 1993; Casotti et al. 1998; Schuchmann 1999; 
GartreU 2000; Gartrell et al. 2000; Nicolson and Fleming 
2003; Downs 2004) and the specialized bill structure of 
some hummingbirds is even considered a result of coevolu- 
tion with the morphology of pollinated flowers (Feinsinger 
and Colwell 1978; Temeles and Kress 2003). 

The adaptation to a new food source like nectar may 
be considered an evolutionary key innovation in the sense 
that it creates an ecological opportunity and promotes 
species proliferation associated with expansion into a pre- 
viously unused niche (Vermel) 1995; Yoder et al. 2010). 
In nectarivorous parrots, this may be particularly true for 
the lories, which were found to be unexpectedly species- 
rich given their age compared to the remaining parrot 
lineages (Schweizer et al. 2011). Ecological opportunity 
may lead to strong directional selection and fast adapta- 
tion (Hunter 1998; Kassen 2009; Yoder et al. 2010). 
Indeed, several morphological specializations of the feed- 
ing tract have been described for nectarivorous parrots, 
which may have been essential for them to effectively feed 
on nectar and poUen. The lories in particular appear to 
have gastrointestinal tracts highly adapted to nectarivory 
(cf Gartrell and Jones 2001). Both they and Lathamus 
have muscular tongues with a brush tip allowing them to 
rapidly harvest nectar (Churchill and Christensen 1970; 
Giintert and ZiswUer 1972; Richardson and WooUer 1990; 
GartreU and Jones 2001). It was further reported that lor- 
ies have shortened intestines, size-reduced gizzards with 
reduced muscularity and koilin layers as well as special 
adaptations in the esophagus, proventriculus and intestine 
(Giintert 1981; Richardson and WooUer 1990). The Lori- 
culus species analyzed so far and Lathamus both shared 
some of the adaptations of the lories in the esophagus, 
proventriculus and intestine (Giintert 1981). However, 
they were found to have comparatively more muscular 
gizzards and longer intestines, probably allowing them to 
feed on hard food like insects or seeds (Giintert and 
Ziswiler 1972; Giintert 1981; GartreU 2000; GartreU et al. 
2000). The morphological and ecological similarities 
among the different nectarivorous parrot groups may 
indicate parallel evolution driven by natural selection. 



M. Schweizer ef al. 

However, data on the morphometries of the digestive 
tract of parrots have never been analyzed with correction 
for phylogenetic non-independence, and the putative 
adaptations to nectarivory as described above have never 
been statistically assessed. AdditionaUy, the comparisons 
of Richardson and WooUer (1990) and GartreU et al. 
(2000) were based on limited taxon sampling. 

In this study, we therefore tested whether the morpho- 
logical variation found in continuous traits of the 
digestive tract of the nectarivorous parrots reflects pheno- 
type-environment correlations as would be expected if 
some of this variation reflects morphological adaptations 
to a nectarivorous diet. We therefore applied a phyloge- 
netic comparative approach using phylogenetic general- 
ized least squares (PGLS) ANCOVA, with diet as a 
covariate in the model to test for phenotype-environment 
correlations. We moreover tested whether a subset of spe- 
cies in a phylogenetic tree shows a trait shift or evolution- 
ary jump at the base of their clade. Since the lories 
apparently show the strongest dependence on nectar 
among aU nectarivorous parrots, we tested if their diet 
specialization was associated with significantly hastened 
morphological evolution at the base of their radiation. All 
analyses were based on 15 continuous characters of the 
digestive tract and a broad taxon sampling of 78 parrot 
species consisting of representatives of all major groups. 
A phylogenetic hypothesis was obtained using three 
nuclear exons and one mitochondrial gene. 

Material and Methods 

Dissection and morphological 
measurements 

Measurements of the digestive tract were taken from 354 
individual parrots (Table SI). The data are from Giintert 
(1981), complemented with 15 additional species (19 
individuals). Body mass was calculated for every species 
as the mean of the fresh dead or frozen and thawed speci- 
mens dissected. AU weights were rounded to the nearest 
0.1 g. For Micropsitta finschi and Loriculus philippensis, 
body mass values were taken directly from the literature 
(Mayr 1931; Rand and Rabor 1960) and an average value 
was calculated combining literature data and the fresh 
weight of other dissected specimens not used in this 
analysis. AU measurements were either taken to the near- 
est mm using dividers for longitudinal measurements 
(length of esophagus, glandular stomach, and intestine) 
or to the nearest 0.1 mm using calipers under a dissecting 
microscope for the other traits. All the digestive organs 
were eventually fixed in buffered formalin (4%) and are 
stored in the vertebrate collection of the Natural History 
Museum Bern. 
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The digestive tract was removed from specimens and 
spread out by cutting through the mesenteria under a 
watery solution of 0.75% NaCl (isotonic for birds). When 
specimens had been preserved in formalin before dissec- 
tion, it was no longer possible to straighten out the intes- 
tine, and the length had to be measured by means of a 
sewing thread, with which all the curvatures of the intesti- 
nal loops could be followed exactly. As parrots lack Brun- 
ner's glands (glandulae duodenales) and caeca (Ziswiler 
and Farner 1972; Giintert 1981), it is not possible to sub- 
divide the intestine into different sections. Length of 
intestine was measured from the pyloric orifice of the giz- 
zard to the rectal widening into the cloaca. 

The esophagus is tripartite in parrots, consisting of a 
pars cervicalis (beginning at the posterior end of the lar- 
ynx), the ingluvies (crop), and a pars thoracica that leads 
into the glandular stomach. Esophageal glands are 
restricted to the caudal area of pars thoracica. Length of 
the esophagus was defined as the distance from the caudal 
rim of the larynx to the border between esophageal and 
gastric glands (Fig. 1). To determine the extension of 
esophagus glands and the transition between the glandular 
part of the glandular stomach and its intermediate zone 
(see below), the digestive tube was cut open longitudinally. 

The proventriculus or glandular stomach contains the 
gastric compound glands. This glandular part is followed 
by an intermediate zone (zona intermedia), lined with 
mucous glands. Total length of the proventriculus was 
measured from the first compound glands visible through 
the wall of the organ to the entrance into the gizzard. 
The caudal measuring point was the cranial groove, situ- 
ated on the pyloric side of the proventricular tube 
(Fig. 2). The extent of the intermediate zone was com- 
puted as total length minus the glandular part (distance 
between the first and the most caudal compound glands). 

The gizzard or muscular stomach has two opposing 
pairs of antagonistic muscles (Figs. 2, 3). Its inner surface 
is lined with a tough koilin membrane (Akester 1986), 
the cuticle of McLeUand (1979), formed by mucosal 
glands. As external dimensions, we measured gizzard 




Figure 1. Longitudinal section tlirougli tine epithelium at the border 
(arrow) between the esophagus (right) and the proventriculus (left) of 
Psittrichas fulgidus. The compound glands (CG) of the proventriculus 
can be clearly distinguished from the mucous glands (MG) of the 
lower part of the esophagus. 
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height (distance between the cranial and caudal groove), 
gizzard depth (minimum distance between tendineal cen- 
ters on the two flat sides of the gizzard), and gizzard 
width (maximum distance at right angles of gizzard 
height) (Fig. 2). Maximum height at main muscles 
(MHM, thick muscle pair) and maximum height at thin 
muscles (MHT) were measured along the maximal exten- 
sion of each muscle pair. Width of caudoventral thin 
muscle (WTM) was measured from the caudal groove to 
its outermost muscle bundles on the opposite side. Width 
of the lumen plus koilin layer (LWiK, including the dis- 
tinctly visible tunica mucosa) was measured along the axis 
of the maximum height at main muscle MHM. Thickness 
of the two thick muscles (MMT) was calculated as the 
difference between MHM and LWiK. Lumen height 
(MLT) was quantified as the maximum distance between 
the opposite walls of the cranial and caudal sac. 

Phylogenetic analyses 

To control for phylogenetic nonindependence in trait val- 
ues, we reconstructed a phylogenetic hypothesis based on 
partial sequences of the three nuclear genes c-mos, RAG-1, 
and Zenk (second exon) and of the mitochondrial gene 
NADH dehydrogenase 2 (ND2) (Table S2). Pitta and Fal- 
co were used as outgroups and the tree was rooted with 
the latter taxon, but both were subsequently pruned from 
the tree before statistical analyses. Sequences of the three 
nuclear genes were taken from Schweizer et al. (2010) or 
newly generated following the laboratory protocol 
described in that study. ND2 sequences were taken from 
GenBank or generated using the primers MetL and ASNH 
for PGR amplification and sequencing from both sides 
(Tavares et al. 2006). The laboratory methods followed 
Schweizer et al. (2010) using the PGR Protocol of Tavares 
et al. (2006) for ND 2 with the annealing temperature set 
to 53°G. The alignment of the sequences was done manu- 
ally with BioEdit 7.0.5.2 (Hall 1999). We checked individ- 
ual sequences and the whole alignment further for quality 
by searching for apparent stop codons after the translation 
of sequences into amino acids. The final alignment was 
4254 bp in length with 603 bp from c-mos, 1461 bp from 
RAG-1, 1149 bp from Zenk, and 1041 bp from ND2. It 
contained one indel of four amino acids for c-mos, one 
indel of three amino acids, and one indel of one amino 
acid for RAG-1, while for Zenk there were four indels of 
one amino acid and one indel of two amino acids. There 
were no ambiguously aligned amino acids. 

Phylogenetic analyses were conducted with maximum 
likelihood (ML) using RAxML 7.0.4 (Stamatakis 2006). 
The program was run on the Web-server with 100 rapid 
bootstrap inferences with all free model parameters esti- 
mated by the software (Stamatakis et al. 2008). We tested 
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Figure 3. Schematic illustration of a transverse 
section through the gizzards along the median 
plane between the tendineal centres of a 
nectarivorous {Vini australis, left) and a 
granivorous parrot {Neophema chrysostoma , 
right). 



different biologically relevant parameter settings for our 
concatenated data set corresponding to separate models 
of nucleotide substitution for genes and/or codon posi- 
tions as estimated by RAxML. Akaike's information crite- 
rion (AIC) was used as a heuristic indicator for the fit of 
the different parameter settings (Akaike 1974), and using 
separate models of nucleotide substitution for the three 
codon positions of each gene separately was found to be 
the best-fitting model. This best-scoring ML tree with 
branch lengths of the best parameter setting was then 
used for further analyses. 

Phylogenetic generalized least squares 
(PGLS) ANCOVA 

PGLS ANCOVA (Grafen 1989; Garland and Ives 2000) 
was used to estimate the relationships between our depen- 



dent and independent variables. This method incorporates 
a matrbc of variance and covariance into the calculation 
of regression parameters based on the pattern of related- 
ness among species to account for the strong correlation 
in the error term (see also below). The variance-covari- 
ance matrix was calculated using the best-scoring ML 
tree. To assess the strength of phylogenetic signal in our 
data, we adjusted our model to include the parameter A, 
which varies between 0 and 1 (Pagel 1997, 1999; Freckl- 
eton et al. 2002). Values of ?. close to 1 imply that traits 
covary as assumed by a Brownian motion model with the 
original tree recovered, while values of 1 close to 0 imply 
that there is almost no phylogenetic signal in the trait 
data, with the phylogenetic tree for the trait having a sin- 
gle polytomy at the basal node (Freckleton et al. 2002; 
Blomberg et al. 2003; Freckleton 2009). A can be inter- 
preted as having one component of the residuals evolving 
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under a Brownian motion model, while another additive 
component has no phylogenetic correlation (Housworth 
et al. 2004; Lavin et al. 2008). Freckleton et al. (2002) 
showed that a is statistically powerful in detecting 
whether the data show a phylogenetic signal, robust to 
incomplete phylogenetic information and that it performs 
better than Grafen's (Grafen 1989, 1992) p transforma- 
tion. We implemented body mass as the independent var- 
iable, because gut measurements of birds are known to be 
allometrically related to body mass (e.g., Ricklefs 1996; 
Lavin et al. 2008), and the 15 morphometric distances 
described above were used as dependent variables. For all 
statistical analyses, both the independent and dependent 
variables were natural-log transformed. 

To test for significant phenotype-environment correla- 
tions between the different traits in the digestive tract and 
nectarivory, we used diet (nectarivory of the lories, 
Brotogeris, Lathamus, and Loriculus versus the more gen- 
eral diets of other parrots) as a covariate in the model. 
We considered increasingly complex models, beginning 
with simple allometry between the dependent and inde- 
pendent variables, then diet was included as covariate 
using ANCOVA with different intercepts but the same 
slope and finally using ANCOVA with different intercepts 
and different slopes (i.e., diet body mass interactions). 
PGLS ANCOVA including the parameter 1 were fitted by 
maximizing the restricted log-likelihood. To check for 
heteroskedasticity, plots of the residuals versus the fitted 
values were investigated. Akaike's information criterion 
(AIC) was used as a heuristic indicator for the fit of the 
different models (Akaike 1974) and we considered an 
increase in model-fit as significant when the reduction in 
AIC score in a more complex model was >4 (Burnham 
and Anderson 2002). We also compared the AIC of the 
models accounting for phylogenetic nonindependence 
with normal general least square approaches, which imply 
full independence of the data. All statistical analyses were 
performed in R 2.9.0 (R Development Core Team 2009) 
using the ape and nlme packages (Paradis et al. 2004; 
Pinheiro et al. 2009). 

Evolutionary trait sHiifts 

We used an additional statistical approach to test for a 
significant shift in trait values at the base of a particular 
clade that cannot be explained by Brownian motion 
(Appendbc SI). The underlying model for character evo- 
lution is a Brownian motion with sudden jumps which 
represent rapid changes. This so-called Levy model has 
been used by many authors including Huelsenbeck et al. 
(2000), Uyeda et al. (2011) and Landis et al. (2012). If 
one assumes the null-hypothesis of no jumps this leads to 
a test statistic whose distribution is known and which 
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thus allowed us to determine P-values in the context of a 
significance test. The basic model is a generalized linear 
model (GLM), similar to that described in Martins and 
Hansen (1997) and in Garland and Ives (2000): 

y = APi2+bP3 + e (1) 

The vector y in equation (1) (see also Appendix SI) 
represents, for each of the N species, the logarithmic val- 
ues of the measured traits which are considered as depen- 
dent variables. The matrix A consists of a column of ones 
(yielding the intercept of the regression line) and a col- 
umn of logarithmic weights. The linear dependence 
between the logarithms of weight and measured trait for 
each species models the assumption of an allometric rela- 
tionship between these quantities. The normally distrib- 
uted error term e reflects the hypothesis that under 
absence of selection the measured trait evolves according 
to a geometric Brownian motion (and thus the logarithm 
of the trait value follows a standard Brownian motion). 
The estimated phylogenetic relationship between the spe- 
cies again induces a strong correlation in the error term. 
In fact, the (i, j)-th element in the variance-covariance 
matrix S, corresponding to two species i and is propor- 
tional to the total length of the shared branches in the 
phylogenetic tree from root to the last common ancestor 
of i and j. The variance-covariance matrix was again cal- 
culated using the best-scoring ML tree. 

The unknown factor of proportionality a represents the 
drift speed of the Brownian movement and has to be esti- 
mated from the data. A cr-value of, for example, two 
would indicate that the relative rate of trait-change under 
neutral evolution is twice as fast as the rate of change in 
the genes used to estimate the variance-covariance matrix 
based on a phylogenetic tree. 

Selection pressure on a subset of K species will result in 
a disproportionate change in the intercept value /? for 
these species, a change that cannot be explained by 
Brownian motion (neutral evolution) alone. The additive 
term in equation (1) containing the parameter /?3 models 
this possibility of selection pressure for the given subset 
of species. The null-hypothesis of no selection can be 
written as /J3 = 0 and the alternative hypothesis corre- 
sponds to /?j / 0. Large values of \ Ps] wiU support the 
alternative hypothesis and result in large absolute values 
of the test statistic T defined in (3). Under the null- 
hypothesis (absence of selection), the distribution of the 
test statistic is a f-distribution, see Theorem 1. Its proof is 
based on the technique of restricted least squares; the 
only nonstandard feature in our situation is the presence 
of correlation, called heteroskedasticity, in the error term. 
In order to reduce our model to the standard situation in 
restricted least squares, we have first to de-couple the 
error terms. The known distribution of the test statistic 
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allowed us to report two-sided P-values. Low P-values 
support a rejection of the null hypothesis. 

We used this approach to test if the lories show a trait 
shift or an evolutionary jump at the beginning of their 
radiation in the 15 measured continuous characters of the 
digestive tract (see above). 

Results 

Phylogenetic relationships 

Our phylogenetic tree from the maximum likelihood 
analyses was in good agreement with Schweizer et al. 
(2010) (Fig. 4). 

The lories were revealed to be a robustly supported 
monophyletic clade clustering as the sister group of Melo- 
psittacus undulatus. The relationships within them were 
also highly supported with the exception of the mono- 



phyly of the genus Trichoglossus and the position of Lorius 
and Chalcopsitta. The division of the lories into two 
clades is in agreement with Wright et al. (2008). The 
position of Loriculus as the sister group of Agapornis was 
also robustly resolved and in agreement with other studies 
(Wright et al. 2008; Schweizer et al. 2010, 2011). Within 
Platycercini, the sister group relationship of Lathamus to 
a clade consisting of Prosopeia, Cyanoramphus and 
Eunymphicus was highly supported in congruence with 
other studies (Schweizer et al. 2010, 2011, 2013; Joseph et al. 
2011). The sister group relationship of Brotogeris within the 
Arini was in congruence with Tavares et al. (2006) and 
Wright et al. (2008), but not supported in our data. 

PGLS ANCOVA 

For linear body dimensions, a scaling with the 0.33 power 
of the body mass can be expected (Schmidt-Nielsen 




Figure 4. Best-scoring maximum likelilnood tree including bootstrap values above 70% indicated at nodes. TInis tree was used for tine formulation 
of the liypothesis of the variance-covariance matrix of the error terms in the regression analyses. Nectarivorous species are shown in red. 
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1984). Indeed, we found scaling factors for all linear 
dimensions to lie around this expected value, ranging 
between 0.18 and 0.46. With the exception of gizzard 
lumen width including the koilin layer (LWiK), all other 
measurements showed a significant phylogenetic signal as 
indicated by lower AIC when phylogenetic information 



was incorporated in the model as compared to a general- 
ized least squares regression (Table 1). 

For eight traits, a PGLS model including diet as a co- 
variate with different intercepts but the same slope was 
considered to be the best-fitting model, while a simple al- 
lometric model best explained the data for the remaining 



Table 1. Akaike's information criterion (AIC) for standard regressions (GLS) and phylogenetic generalized least squares (PGLS) including the fitted 
1-values for the different models considered for each trait of the digestive tract. Increasingly complex models were tested, beginning with simple 
allometry between the dependent and independent variables, then diet (nectarivory) was included as covariate using ANCOVA with different inter- 
cepts but the same slope (Mass + Food) and finally using ANCOVA with different intercepts and different slopes (Mass x Food). 
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ividii 1 ruuu 


— 52 258 


—58 246 


0 652 




N/l^cc y FnnH 
ividjj A ruuu 


—47 040 


—52 899 


0 661 


C-i\~7~77i'rr\ Honth 
\j\z.i.a\\i utrULii 


'-[imnio ^llnmo+R/ 
-jiiiiuitr diiuiiitriiy 


40 475 


—21 107 


0 914 




-1- FnnH 
ividjj 1 ruuu 


—27 41 3 


—34 577 


0 804 




^l'l^cc V FnnH 

IVIajj A 1 UUU 


—25.41 7 


—31 .727 


0.769 


I\yi3Yimiim niTvariH hoinht main mi it:/"lcic /f\/ll— l^yl^ 
iviciAii 1 lUM 1 yizLzLdiu iitriyiiL cil [iidiii iiiuiiuitr^ \ivinivi/ 


^Nimnlo 3llnmotr\/ 
_)iiii[jitr diiuiiitrLiy 


1 0 205 


—28 377 


0 763 




Mass + Food 


—33.2 1 2 


—36.807 


0.589 




Mass X Food 


-27.998 


-31.641 


0.593 


Gizzard thickness at main muscles (MMT) 


Simple allometry 


169.006 


102.775 


0.938 




Mass + Food 


103.458 


89.283 


0.898 




Mass X Food 


106.376 


91.087 


0.917 


Gizzard lumen width including koilin layer (LWiK) 


Simple allometry 


0.110 


0.104 


0.334 




Mass + Food 


5.594 


4.28 


0.395 




Mass X Food 


10.084 


8.877 


0.398 


Gizzard width at caudoventral thin muscle (WTM) 


Simple allometry 


-37.786 


-47.244 


0.680 




Mass + Food 


-33.593 


-41.924 


0.635 




Mass X Food 


-32.005 


-37.503 


0.667 


Maximum gizzard height at thin muscle (MHT) 


Simple allometry 


-37.340 


-89.566 


0.846 




Mass + Food 


-82.851 


-93.727 


0.759 




Mass X Food 


-76.966 


-88.706 


0.747 


Maximum gizzard lumen at thin muscle (MLT) 


Simple allometry 


-30.160 


-69.159 


0.791 




Mass + Food 


-62.319 


-70.648 


0.699 




Mass X Food 


-56.805 


-65.538 


0.697 


Gizzard mass 


Simple allometry 


130.144 


75.418 


0.880 




Mass + Food 


75.552 


66.876 


0.787 




Mass X Food 


79.187 


70.643 


0.781 
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Trait 


Model 




Value 


SE 


f-value 


P-value 


Length of intestine 


Simple allometry 


Intercept 


4.251 


0.251 


16.957 


0.000 






Slope 


0.393 


0.033 


11.772 


0.000 


Length of esophagus 


Simple allometry 


Intercept 


2.649 


0.083 


31.743 


0.000 






Slope 


0.356 


0.014 


25.71 1 


0.000 


Extension of esophagus glands 


Mass + Food 


Intercept 


0.002 


0.313 


0.006 


0.995 






Slope 


0.321 


0.044 


7.227 


0.000 






Intercept 


0.627 


0.141 


4.431 


0.000 


Length of intermediate zone 


Mass + Food 


Intercept 


0.258 


0.256 


1.008 


0.317 






Slope 


0.456 


0.037 


12.482 


0.000 






Intercept 


-0.221 


0.115 


-4.149 


0.000 


Length of proventriculus 


Simple allometry 


Intercept 


1.404 


0.120 


11.717 


0.000 






Slope 


0.389 


0.019 


20.414 


0.000 


Gizzard height 


Mass + Food 


Intercept 


0.984 


0.116 


8.481 


0.000 






Slope 


0.322 


0.018 


17.775 


0.000 






Intercept 


1.199 


0.053 


4.000 


0.000 


Gizzard width 


Simple allometry 


Intercept 


1.155 


0.158 


7.302 


0.000 






Slope 


0.336 


0.024 


14.166 


0.000 


Gizzard depth 


Mass + Food 


Intercept 


0.434 


0.188 


2.305 


0.024 






Slope 


0.303 


0.027 


11.301 


0.000 






Intercept 


0.864 


0.085 


5.071 


0.000 


Maximum gizzard height at main muscles (MHM) 


Mass + Food 


Intercept 


1.115 


0.159 


7.027 


0.000 






Slope 


0.298 


0.025 


1 1 .988 


0.000 






Intercept 


1.433 


0.074 


4.310 


0.000 


Gizzard thickness at main muscles (MMT) 


Mass + Food 


Intercept 


-0.236 


0.473 


-0.499 


0.619 






Slope 


0.182 


0.063 


2.871 


0.005 






Intercept 


0.713 


0.21 1 


4.501 


0.000 


Gizzard lumen width including koilin layer (LWiK) 


Simple allometry 


Intercept 


0.606 


0.170 


3.568 


0.001 






Slope 


0.354 


0.030 


11.831 


0.000 


Gizzard width at caudoventral thin muscle (WTM) 


Simple allometry 


Intercept 


0.447 


0.155 


2.885 


0.005 






Slope 


0.370 


0.024 


15.300 


0.000 


Maximum gizzard height at thin muscle (MHT) 


Mass + Food 


Intercept 


1.195 


0.122 


9.776 


0.000 






Slope 


0.320 


0.018 


17.920 


0.000 






Intercept 


1.391 


0.055 


3.542 


0.000 


Maximum gizzard lumen at thin muscle (MLT) 


Simple allometry 


Intercept 


1.061 


0.145 


7.293 


0.000 






Slope 


0.348 


0.022 


16.081 


0.000 


Gizzard mass 


Mass + Food 


Intercept 


-4.505 


0.365 


-12.352 


0.000 






Slope 


0.870 


0.053 


16.559 


0.000 






Intercept 


-3.854 


0.165 


3.947 


0.000 



traits (Tables 1, 2, Figs. 5, 6). The specialization to nec- 
tarivory led to a decrease in the extension of the esopha- 
gus glands, though the trait-value of Lathamus was 
seemingly more similar to the non-nectarivorous parrots 
and there was some variation within the lories. Further- 
more, the length of the intermediate zone was found to 
be prolonged in nectarivorous parrots, though this was 
apparently not the case for Brotogeris. 

The nectarivorous species clearly differed from the 
remaining parrots in the measurements of the gizzard, 
with the exception of the following measures: width of 
whole gizzard and its lumen, width at the caudoventral 
thin muscle (WTM) and the maximum lumen at thin 
muscles (MLT). Brotogeris and Lathamus were seemingly 



more like the nonnectarivorous parrots for all gizzard 
traits. This was apparently also the case for Loriculus galg- 
ulus in gizzard mass, gizzard height, gizzard thickness at 
main muscles (MMT), and maximum gizzard height at 
main muscle (MHM). The gizzard measurements of Psit- 
trichas fulgidus showed a clear tendency to be more clo- 
sely associated with the nectarivorous species than with 
the nonnectarivorous parrots. 

In contrast, including diet as covariate in the model 
did not improve model-fit for the length of the esophagus 
or the intestine. While an intermediate value of I was 
revealed for the esophagus length, the length of the intes- 
tine showed a strong phylogenetic signal with a value of X 
close to one (Table 1). When diet was included in the 



2874 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



M. Schweizer ef al. 



Parallel Adaptations to Nectarivory In Parrots 
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In weight 

Figure 5. Natural-log transformed values of the different independent variables against natural-log transformed body masses including the 
regression lines of the best-fitting model. For all these traits, the data were best explained by an allometric relationship between the dependent 
and independent variables. 
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Figure 6. Natural-log transformed values of the different independent variables against natural-log transformed body masses including the 
regression lines of the best-fitting models. For all these traits, a model including the nectarivorous diet as a covariate with different intercepts but 
the same slope was considered as the best-fitting model. 



model for the length of the proventriculus, the AIC was Table 1). However, the more complex model was not 
slightly lower for an ANCOVA with the same slopes but substantially supported. Hence, the length of the proven- 
different intercepts compared to simple aUometry (Fig. 5, triculus is best explained by simple aUometry. 
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Table 3. P-values and drift speed of the test for rate shifts in mor- 
phological evolution at the base of the lory radiation. 



Trait 


P-value 


(T-value 
(drift speed) 


Length of intestine 


n.s. 


0.9667 


Length of esophagus 


n.s. 


0.5011 


Extension of esophagus glands 


2.26E-04 


1.2912 


Length of intermediate zone 


1 .30E-04 


1.0891 


Length of proventriculus 


n.s. 


0.6343 


Gizzard height 


0.0058 


0.6319 


Gizzard width 


0.021 


0.756 


Gizzard depth 


2.99E-04 


0.8304 


Maximum gizzard height at 


0.0254 


0.9091 


main muscles (MHM) 






Gizzard thickness at main 


2.00E-06 


1.7508 


muscles (MMT) 






C-i\77^rr\ liimpn wirlth inrliirlinn 

VJ IZ-Z-OIU ILIIIICril VVIULII lllv.lLJLJIII^J 




1 .2635 


koilin layer (LWiK) 






Gizzard width at caudoventral thin 


n.s. 


0.8016 


muscle (WTM) 






Maximum gizzard height at 


0.0011 


0.5496 


thin muscle (MHT) 






Maximum gizzard lumen at 


0.0091 


0.6638 


thin muscle (MLT) 






Gizzard mass 


6.86E-05 


1.5585 



Evolutionary trait shifts 

For ten traits, we found a significant shift; in trait evolu- 
tion at the base of the lory radiation (Table 3). These 
included all traits for which a model with diet as covari- 
ate was the best-fitting model in the PGLS ANCOVA's. In 
addition to these traits, a significant shift in trait evolu- 
tion at the base of the lories was found for gizzard width 
and MLT. When correcting for multiple comparisons, 
however, the shift in trait evolution was no longer signifi- 
cant for MLT and neither for gizzard height, gizzard 
width and MHM. 

Discussion 

Phenotype-environment correlations 

The gastrointestinal tract may be considered one of the 
major interfaces between an organism and its environment, 
mediating their interactions (Karasov 1990). Several inter- 
specific studies on other bird groups have shown that the 
size of the gut is related to diet and that the morphology of 
the gastrointestinal tract often reflects the birds' feeding 
strategies (Ricklefs 1996; Battley and Piersma 2005; Cavie- 
des-Vidal et al. 2007; Lavin and Karasov 2008). Moreover, 
a direct influence of feeding strategies on the structures, 
functionality and physiology of the digestive tract has been 
shown in other vertebrates such as mammals (e.g., Schieck 



and Millar 1985; Korn 1992; Lovegrove 2010), amphibians 
and reptiles (e.g., Stevens and Hume 1995; O'Grady et al. 
2005), and fish (e.g., German and Horn 2006; Wagner et al. 
2009). However, it has to be considered that biological 
structures are not only fine tuned to their functional 
demands by natural selection but are also influenced by 
phylogenetic history and biochemical and mechanical con- 
straints (Raia et al 2010). 

Within parrots, we have shown here that nectarivory is 
associated with reduced extension of the glands in the 
lower part of the esophagus below the crop (Pars thoraci- 
ca). These glands produce a mucous secretion which 
helps hard ingesta to glide through the glandular stomach 
(proventriculus) and reduces the risk of mechanical dam- 
age to the latter (Giintert 1981). It can be expected that 
parrot species eating exclusively soft or liquid food evolve 
reductions of these glands within the esophagus, and our 
data corroborate this. On the other hand, we could not 
find any indication that the nectarivorous parrots have a 
longer esophagus as was proposed for the lories by 
Giintert (1981). 

Between the glandular stomach (proventriculus) and 
the gizzard (muscular stomach, ventriculus), there is an 
intermediate zone characterized by the absence of com- 
pound glands of the former and absence of the koilin 
layer of the latter (Ziswiler and Farner 1972). This inter- 
mediate zone has the function of a storage space, where 
the proteolytic enzyme pepsin from the proventriculus 
can react with ingesta (Giintert 1981). We found the nec- 
tarivorous parrots to have a longer intermediate zone 
compared to the remaining parrots, even though the trait 
values of one nectarivorous taxon, Brotogeris, were more 
similar to the non-nectarivorous species. The intermediate 
zone appears to play an important role in the digestion of 
poUen (Giintert 1981) and a prolonged intermediate zone 
may be an adaptation to optimize the extraction of amino 
acids from pollen grains. Pollen grains have a high pro- 
tein content, with their interior (protoplast) consisting of 
diverse amino acids (van Tets and Hulbert 1999; Gartrell 
and Jones 2001). Acidifications of pollen grains in the 
proventriculus may be important so that their contents 
can be extruded and digested, whereas mechanical break- 
up of pollen grains in the gizzard does not seem to be 
important (Gartrell and Jones 2001). The amount of 
energy extracted from meals can be enhanced by increas- 
ing the retention time (McWhorter et al. 2009). A pro- 
longation of the intermediate zone can thus increase the 
rate of protein digestion as it increases the retention time 
of pollen grains, which seems to be the case in nectarivor- 
ous parrots. PoUen ingestion may require less energy than 
feeding on insects as an additional amino acid supply, 
because nectarivorous birds will encounter pollen while 
feeding on nectar (Nicolson and Fleming 2003). 
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Besides being a storage organ, the gizzard functions as 
an organ of mechanical digestion, the site of preliminary 
acid proteolytic digestion, and a filter for indigestible mate- 
rial (Ziswiler and Farner 1972). During contraction of the 
gizzard, the thick muscles close up, narrowing the lumen 
to a thin cleft and forcing the contents into two pouches 
(cranial and caudal sac, cf Fig. 2) that lie under the thin 
muscles (McLelland 1979). Species feeding on a soft diet 
do not need the grinding function to break down their 
food, and can be expected to evolve reduced gizzard mus- 
culature (Steinbacher 1934; McLelland 1979). Indeed, the 
nectarivorous parrots differed from the remaining parrots 
by having less developed gizzard muscles. In contrast, a 
simple aUometric relationship best explained the width of 
the whole gizzard and its lumen as well as its width at 
the caudoventral thin muscle and the maximum lumen at 
the thin muscles. The thin muscles act antagonistically to 
the main muscles and have no grinding function. There- 
fore, they are not expected to be developed more strongly 
in species relying on the grinding function of the gizzard. 
In congruence with our results, Richardson and WooUer 
(1990) also found two species of lories [Glossopsitta por- 
phyrocephala, Trichoglossus haematodus) to have smaller 
and less muscular gizzards than four nonnectarivorous 
parrot species of Australia {Melopsittacus undulatus, Bar- 
nardius zonarius, Neopsephotus bourkii, Platycercus ictero- 
tis). Interestingly, we found that the reportedly mainly 
frugivorous Pesquet's parrot Psittrichas fulgidus (Collar 
1997) shared a similarly reduced muscularity with the nec- 
tarivorous parrots. On the other hand, the blue-crowned 
hanging-parrot Loriculus galgulus did not show an overall 
reduced muscularity. This may correspond to the higher 
amount of seeds in its diet compared to L. philippensis 
(Homberger 1980). Similarly, gizzard measurements of 
Lathamus discolor and Brotogeris jugularis clustered with 
the non-nectarivorous parrots, and these two species were 
not found to have overall reduced gizzard muscularity 
either. Other studies also found Lathamus discolor to have 
retained the muscular gizzard of a granivorous species 
(Giintert and Ziswiler 1972; Gartrell et al. 2000). This may 
allow this species to feed on harder food when nectar and 
pollen are rare (Gartrell et al. 2000). Gizzard dimensions 
also vary in other passerines according to diet, with longer 
gizzards in seed- than in fruit- and insect-eaters and 
thicker muscular and glandular layers in insect- compared 
to fruit- and seed-eaters (Ricklefs 1996). Smaller gizzards 
with a reduced muscularity were also found in the necta- 
rivorous honeyeaters (Meliphagidae) compared to similar- 
sized passerines (Richardson and WooUer 1986); however, 
phylogenetic nonindependence was not controlled for in 
that study. 

Chemical digestion of food principally takes place in 
the intestine (Ziswiler and Farner 1972). Richardson and 



WooUer (1990) found two species of lories (Glossopsitta 
porphyrocephala, Trichoglossus haematodus) to have 
shorter intestines compared with four non-nectarivorous 
parrot species of Australia [Melopsittacus undulatus, Bar- 
nardius zonarius, Neopsephotus bourkii, Platycercus ictero- 
tis). This was explained as a consequence of sugars in 
nectar needing less processing in the intestine than other 
food. In contrast, we could not find any indication for 
shorter intestines of nectarivorous parrots. In a broad 
comparative study of birds, Lavin et al. (2008) did not 
find any significant effect of diet on small intestine 
length either (Lavin et al. 2008). However, this result is 
only partly comparable with ours, since we measured 
the whole lengths of the intestine owing to the difficulty 
of clearly distinguishing the small and the large intestine 
in parrots due to lack of caeca (Ziswiler and Farner 
1972; Giintert 1981). In addition to lengths, intestine 
function depends inter alia on volume, surface area, 
viUi, and microviUi area as well as enzymatic activity 
(Ricklefs 1996; Lavin et al. 2008; McWhorter et al. 
2009). Moreover, the efficiency of digestion in the intes- 
tine may be influenced by the passive absorption of hy- 
drosoluble compounds through the paraceUular pathway. 
This is prominent in birds and may be especially impor- 
tant for nectarivores because they have to deal with 
large amounts of sugar in their diet (Karasov and Cork 
1994; Napier et al. 2008; McWhorter et al. 2009). In 
general, birds have a lower nominal surface area of the 
intestine and a shorter smaU intestine as weU as shorter 
digestive retention times than mammals; however, their 
higher passive absorption compared to mammals may 
compensate for this (McWhorter et al. 2009). This may 
render predictions about the intestine dimensions in 
relation to diet more difficult. 

In conclusion, our analyses showed that nectarivorous 
parrots differ, after correction for phylogenetic noninde- 
pendence, from the remaining parrots in several traits 
of the digestive tract. Hence, we uncovered significant 
phenotype-envirormient correlations for the prolonga- 
tion of the intermediate zone, the reduction of gizzard 
muscularity and the reduction of glands in the esopha- 
gus. The similarity in these trait features among some 
of the different nectarivorous groups is an indication of 
parallel evolution under the same or similar environ- 
mental conditions, that is, the shift to a nectarivorous 
diet, and implies that natural selection was the main 
driving force (cf Losos et al. 1998; Schluter et al. 2004; 
Colosimo et al. 2005). Moreover, functional consider- 
ations suggest that the adaptations in the intermediate 
zone of nectarivorous parrots (probably except Brotoge- 
ris) allow them to rely effectively on nectar as a food 
source, and thus implying evidence for trait utility (Sch- 
luter 2000). 
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Phenotypic flexibility of the gastrointestinal 
tract 

The size, structure, and functional characteristics of the 
gastrointestinal tract of birds can reversibly change 
within the lifetime of a bird (phenotypic flexibility, 
sensu Piersma and Drent (2003)) as a fast adaptive 
response to current functionality demands caused by 
environmental changes or circannual endogenous control 
(Starck 1999a,b; Piersma and Drent 2003; Starck and 
Rahmaan 2003; Battley and Piersma 2005; McWhorter 
et al. 2009). As pointed out by Lavin et al. (2008), com- 
parative studies like ours have the limitation that species 
were not analyzed under common-garden conditions 
and thus it is not possible to assess to what extent the 
variation found among species is influenced by pheno- 
typic flexibility and plasticity at the individual level. 
However, the individuals analyzed in this study all stem 
from captivity, where more stable conditions than in 
nature can be expected, thus mirroring a common-gar- 
den experiment. Furthermore, the inclusion of several 
individuals for most species and the wide range of body 
sizes among species considered certainly minimized the 
effect of intraspecific variation. There is additionally 
some evidence that phenotypic flexibility of the gastroin- 
testinal tract is limited in parrots (cf Giintert 1981). All 
individuals of Lathamus for example analyzed in this 
study were fed with a nectar-alternative, but they 
retained the partly muscular gizzard similar to that of a 
granivorous species, and the features of their gastrointes- 
tinal tract did not appear to differ from wild specimens 
analyzed by Gartrell et al. (2000). Nevertheless, further 
studies, preferably on wild birds, are needed to docu- 
ment the interplay between natural selection, plasticity 
and potential flexibility in features of the digestive tract 
in different parrot species. 

Evolutionary trait shifts and species 
proliferation 

We found that in the lories, the diet shift to nectarivory 
was associated with a significant shift in morphological 
evolution, chiefly of several gizzard traits, at the base of 
their radiation, implying a trait shift or evolutionary jump. 

The lories have diversified into an exceptionally spe- 
cies-rich clade (Schweizer et al. 2011) and their diet shift 
might thus be considered to be an evolutionary key inno- 
vation. Nectar may have provided a spatially widespread 
underutilized niche and this may have allowed the lories 
to expand their ranges and to colonize even remote oce- 
anic islands, which may have fostered allopatric specia- 
tion. Even today, congeneric species of the lories generally 
do not overlap geographically (Collar 1997). Sympatry 



within genera is found in eastern Australia, New Guinea 
and Wallacea, regions with a complex and composite 
environmental and geological history with several poten- 
tial vicariance opportunities in the past (Hall 2002; Essels- 
tyn et al. 2009; Byrne et al. 2011; Deiner et al. 2011). 
However, this ecological expansion was not followed by 
further significant ecological specializations within the 
radiation of the lories. Similar to honeyeaters, the other 
highly nectarivorous and species-rich bird-group of Aus- 
tralasia (Newton 2003), lories are generalized flower visi- 
tors and their ecological relationships with plants are not 
as specialized as those of hummingbirds or sunbirds 
(Fleming and Muchhala 2008). Avian pollinator assem- 
blages differ regionally and the evolutionary specializa- 
tions between nectar-feeding birds and their food-plants 
are strongest in the Neotropics, decreasing through Africa 
and South Asia to Southeast Asia and Australasia (Flem- 
ing and Muchhala 2008). The co-evolution of specialized 
plant-pollinator relationships may take time. While the 
hummingbirds split from their closest relatives in the 
Eocene (about 50 Ma) or even earlier (Brown et al. 2008; 
Pratt et al. 2009) with a major radiation after 20 Ma (Jetz 
et al. 2012), the lories split from their common ancestor 
with Melopsittacus only in the middle Miocene (about 
15 Ma, Schweizer et al. 2011). Hummingbirds certainly 
had more time to co-evolve with plants than did lories. 
However, the similarly specialized sunbirds are likely to 
be younger than hummingbirds, and the evolutionary 
diversification of the generalist honeyeaters started at a 
similar age (Eocene) (Barker et al. 2004). Thus, explana- 
tions other than time may account for the low ecological 
specialization of Australasian nectarivorous birds com- 
pared with sunbirds and hummingbirds. Specific interac- 
tions between plants and pollinators are only likely to 
evolve when floral resources are spatially and temporally 
predictable (e.g., Waser et al. 1996). While this seems to 
be the case in the Neotropics, flowering of eucalypts in 
Australia varies in space and time, and trees in lowland 
and montane Papua New Guinea commonly have non- 
annual flowering patterns (Fleming and Muchhala 2008). 
Birds feeding on them may hence not be able to afford to 
specialize. This may account for the low specialization in 
plant-pollinator relationships of Australasian nectarivor- 
ous birds and may explain the lack of evolution of plant- 
specific specializations among the lories after their shift to 
a nectarivorous diet. 

In conclusion, the key innovation of the lories allowed 
an expansion into a new adaptive zone and we hypothe- 
size that the subsequent species proliferation may have 
essentially been nonadaptive through allopatric speciation. 
The lories may thus be considered an example of a non- 
adaptive radiation (Rundell and Price 2009). It is possible 
that the ecological opportunity provided by their key 
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innovation did not trigger an adaptive radiation because 
of the unpredictable nature of the new resource. The key 
innovation nevertheless promoted significant lineage 
diversification through allopatric partitioning of the same 
broad new niche. Although other parrot groups switched 
to a nectarivorous diet, this did not increase their diversi- 
fication rates and species richness compared to other par- 
rots (Schweizer et al. 2011). Various factors may have 
inhibited an increased rate of cladogenesis in Lathamus, 
Brotogeris, and Loriculus following their change to a nec- 
tarivorous diet. Such factors can include developmental 
or genetic constraints, but also ecological circumstances 
like interspecific competition or the lack of opportunities 
for allopatric speciation. Hence, an evolutionary innova- 
tion does not necessarily lead to increased diversification 
(Vermeij 2001; Price et al. 2010). The question as to 
which factors hampered increased species proliferation in 
other nectarivorous parrots will be an interesting avenue 
for future research. 
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